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Abstract 

We study the complexity of solving the generalized MinRank problem, i.e. computing the 
set of points where the evaluation of a polynomial matrix has rank at most r. A natural algebraic 
^ representation of this problem gives rise to a determinantal ideal: the ideal generated by all minors 

O of size r + 1 of the matrix. We give new complexity bounds for solving this problem using Grobner 

bases algorithms under genericity assumptions on the input matrix. In particular, these complexity 
,—1 bounds allow us to identify families of generalized MinRank problems for which the arithmetic 

J> complexity of the solving process is polynomial in the number of solutions. We also provide 

an algorithm to compute a rational parametrization of the variety of a 0-dimensional and radical 
system of bi-degree (D, 1). We show that its complexity can be bounded by using the complexity 
bounds for the generalized MinRank problem. 

1 Introduction 

We focus in this paper on a problem which admits a natural algebraic formulation as a structured 
• "-j system of polynomial equations: 

^ Generalized MinRank Problem: given a field IC, a « x m matiix whose entries are polynomials 

^ of degree D over lC[.xi, . . . ^Xk\, and r < min{n,m) an integer, compute the set of points at which the 

evaluation of the matrix has rank at most r. 

Being able to estimate precisely the complexity of this problem (which is known to be NP- 
complete when IC is a finite field (5]) is of first importance for applications. In Cryptology, the 
security of several multivariate cryptosystems relies on the difficulty of solving the classical MinRank 
problem (i.e. when the entries of the matrix are linear ETlfmi ). In coding theory, rank-metric codes 
can be decoded by computing the set of points where a polynomial matrix has rank less than a given 
value |[23l [TTll . Also, in Geometry and Optimization the critical points of a map are defined by the 
rank defect of its Jacobian matrix. Moreover, this problem is also underlying other problems from 
Symbolic Computation (for instance solving multi-homogeneous systems, see e.g. [13|). 

To study the Generalized MinRank problem, we consider the algebraic system of all (r+ 1)- 
minors of the input matrix. Indeed, these minors simultaneously vanish on the locus of rank defect 
and hence give rise to a section of a determinantal ideal. 
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Several solving tools can be used to solve this algebraic system by taking profit of the underlying 
structure. For instance, geometric resolution flTl can use the fact that these systems can be evalu- 
ated efficiently. Also, recent works on homotopy methods show that numerical algorithms can solve 
determinantal problems |[25l . 

In this paper, we focus on Grobner bases algorithms. A representation of the locus of rank defect 
is obtained by computing a lexicographical Grobner basis by using the algorithms F5 ||9l and FGLM 
ifTOll . Indeed, experiments suggest that these algorithms take profit of the determinantal structure. 
The aim of this work is to give an explanation of this behavior from the viewpoint of asymptotic 
complexity analysis. 



Related works 



An important related theoretical issue is to understand the algebraic structure of the ideal C IK[?7] 
(where U is the set of variables {mi 1 , . . . , m„,„j}) generated by the (r + l)-minors of the matrix: 



/ "1,1 



The ideal has been extensively studied during last decades. In particular, explicit formulas for 
its degree and for its Hilbert series are known (see e.g. |[T6l Example 14.4.14] and 0), as well as 
structural properties such as Cohen-Macaulayness and primality |fT9ll20l . 

In Cryptology, Kipnis and Shamir have proposed a multi-homogeneous algebraic modeling which 
can be seen as a generalization of the Lagrange multipliers [21] and is designed as follows: a poly- 
nomial nxm matrix ^ G ]K[X]"^'" (where X denotes the set of variables {x\ ,xi(}) has rank less 
than r + 1 if and only if the dimension of its right kernel is greater than m — r — 1. Consequently, by 
introducing r(m — r) fresh variables 1 , . . . ,yr,m-r, we can consider the system of bi-degree (D, 1) in 
K[xi,... ^Xk-,y 1,1, ■ ■ ■ ,yr,m-r] defined by 
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If {x\, . . . ,Xk,yi,\, ■ ■ ■ ,yr,m-r) is a solution of that system, then the evaluation of the matrix ^ at 
the point (jci , . . . ,Xk) has rank less than r + 1 . 

In 1 12 1, the case of square linear matrices is studied by performing a complexity analysis of the 
Grobner bases computations. In particular, this investigation showed that the overall complexity is 
polynomial in the size of the matrix when the rank defect n — ris constant. This theoretical analysis is 
supported by experimental results. The proofs were complete when the system has positive dimension, 
but depended on a variant of a conjecture by Froberg in the 0-dimensional case. 
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Main results 



We generalize in several ways the results from |[T2]| . where only the case of square linear matrices 
was investigated. The goal of this paper is to obtain complexity bounds for Grobner bases algorithms 
when the input system is the set of (r+ l)-minors of an xm matrix whose entries are polyno- 
mials of degree D with generic coefficients. By generic, we mean that there exists a non-identically 
null multivariate polynomial h such that the complexity results holds when this polynomial does not 
vanish on the coefficients of the polynomials in the matrix. Therefore, from a practical viewpoint, the 
complexity bounds can be used for applications where the base field IC is large enough: in that case, 
the probability that the coefficients of do not belong to the zero set of h is close to 1 . 

We start by studying the homogeneous generalized MinRank problem (i.e. when the entries of 
^ are homogeneous polynomials) and by proving an explicit formula for the Hilbert series of the 
ideal generated by the (r+ 1) -minors of the matrix The general framework of the proofs is 
the following: we consider the ideal C K[U] generated by the (r+ l)-minors of a matrix ^ = 
(ujj) whose entries are variables. Then we consider the ideal &r = + {gi, ■ ■ ■ ,gnm) C ]K[?7,X], 
where the polynomials gi are quasi-homogeneous forms that are the sum of a linear form in K[U] 
and of a homogeneous polynomial of degree D in K[X]. If some conditions on the gi are verified, by 
performing a linear combination of the generators there exists /i.i, . . . ,/„,„, G K[X] such that 



Then we use the fact that (^^ + ("1.1 — /i.i, • • • )"ji,m — /«,m)) nlC[X] = J^r to prove that properties 
of generic quasi-homogeneous sections of transfer to J^^ when the entries of the matrix ^ are 
generic. This allows us to use results known about the ideal to study the algebraic structure of J^^. 

We study separately three different cases: 

• k> {n — r){m — r). Under genericity assumptions on the input, the solutions of the generalized 
MinRank problem are an algebraic variety of positive dimension. The complexity results were 
only proven for D = 1 and « = m in [12J. We generalize here for any D G N. 

• k = {n — r){m — r). This is the — dimensional case, where the problem has finitely-many 
solutions under genericity assumptions. The results in [12J were only stated for D = 1 and 
n = m, and they depended on a variant of Froberg's conjecture. In this paper, we give complete 
proofs for D G N which do not rely on any conjecture. 

• k < {n — r){m — r). In the over-determined case, we still need to assume a variant of Froberg's 
conjecture to generalize the results in |[T2l . 

In particular, we prove that, for ^ > {n — r){m — r), the Hilbert series of the coordinate ring of 
is the power series expansion of the rational function 




where Ar{t) is the r xr matrix whose (/,7)-entry is Y,k ('\ {" k^)^'^- Assuming w.l.o.g. that m <n, 
we also prove that the degree of is equal to 



m—r— 1 

DEG(^,) =d("-'")('"-'") Y{ 

i=0 
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These explicit formulas permit to derive complexity bounds on the complexity of the problem. 
Indeed, one way to get a representation of the solutions of the problem in the 0-dimensional case 
is to compute a lexicographical Grobner basis of the ideal generated by the polynomials. This can 
be achieved by using first the ^5 algorithm [9] to compute a Grobner basis for the so-called grevlex 
ordering and then use the FGLM algorithm lITOl to convert it into a lexicographical Grobner basis. 
The complexities of these algorithms are governed by the degree of regularity and by the degree of 
the ideal. 

Therefore the theoretical results on the structure of J^^ yield bounds on the complexity of solv- 
ing the generalized MinRank problem with Grobner bases algorithms. More specifically, when k = 
[n — r){m — r) and under genericity assumptions on the input polynomial matrix, we prove that the 
arithmetic complexity for computing a lexicographical Grobner basis of is upper bounded by 



Cl) 



n\( m \/Dreg+^\ ,wnl=r^^^^3 



where 2 < ft) < 3 is a feasible exponent for the matrix multiplication, and 

D^gg = Dr(m - r) + (D - l)k+\. 

This complexity bound permits to identify families of Generalized MinRank problems for which 
the number of arithmetic operations during the Grobner basis computations is polynomial in the num- 
ber of solutions. 

In the over-determined case (i.e. k < {n — r){m — r)), we obtain similar complexity results, by 
assuming a variant of Froberg's conjecture which is supported by experiments. 

Finally, we show that complexity bounds for solving systems of bi-degree (D, 1) can be obtained 
from these results on the generalized MinRank problem. We give an algorithm whose arithmetic 
complexity is upper bounded by 

H(::"')r'::'")"- 

for solving systems of rix + riy equations of bi-degree {D,l) in lC[xi, . . . ,x„^,ji, . . . ,j„^] which are 
radical and 0-dimensional. 




Organization of the paper 

Section[2]provides notations used throughout this paper and preliminary results. In Section[3} we show 
how properties of the ideal generated by the (r + 1) -minors of ^ transfer to the ideal Then, the 
case when the homogeneous Generalized MinRank Problem has non-trivial solutions (under gener- 
icity assumptions) is studied in Section [4] Section [5] is devoted to the study of the over-determined 
MinRank Problem (i.e. when k < {n — r){m — r)). Then, the complexity analysis is performed in Sec- 
tion [6] Some consequences of this complexity analysis are drawn in Section [7] Experimental results 



are given in Section 7.4 and applications to the complexity of solving bi-homogeneous systems of 



bi-degree (D, 1) are investigated in Sections] 



2 Notations and preliminaries 

Let K be a field and K be its algebraic closure. In the sequel, n, m, r and k and D are positive integers 
with r < m < n. For S N, Mon{d,k) denotes the set of monomials of degree d in the polynomial 
ring K[xi ,Xk]. Its cardinahty is #Mon{d,k) = ('^"]^'') • 



4 



We denote by o the set of parameters {am^^ '■ ^ < i < n,l < j < m,t ^ Mon{D,k)}. The set of 
variables {ujj : I < i < n, I < j < m} (resp. {xi , . . . ,Xk}) is denoted by U (resp. X). 
For 1 < / < ?i, 1 < 7 < m, we denote by fi j G lC(a)[X] a generic form of degree D 



teMon{DM) 

Let C ]K(o)[X] be the ideal generated by the (r + 1) -minors of the nxm matrix 

• • • fn 

and C K{a)[U,X] be the determinantal ideal generated by the (r+ l)-minors of the matrix 

/ui I ... Ui „ 



We define Jy as the ideal 'Sy + - fi,j)\<i<n.\<j<m C IK(a) [?7,X]. Notice that Jy = {uij - 
Aj)i<i<nA<j<m C K{a)[U,X]. Therefore, = JynK{a)[X]. 

By slight abuse of notation, if / is a proper homogeneous ideal of a polynomial ring K[X], we call 
Hubert series of / and we note HS/ G Z[[?]] the Hilbert series of its coordinate ring with the 

grading defined by deg(;t,) = 1 for all /: 

HS/(0= £dimK(lC[X],///rf)f^ 

d>0 

where denotes the vector space of homogeneous polynomials of degree d and Id = I nK[X]d. 
We call dimension of / the KruU dimension of the quotient ring 



Quasi-homogeneous polynomials. 

We need to balance the degrees of the entries of the matrix with the degrees of the entries of 
This can be achieved by putting a weight on the variables Ujj, giving rise to quasi-homogeneous 
polynomials. A polynomial / G K[?7,X] is called quasi-homogeneous (of type (D, 1)) if the following 
condition holds (see e.g. [18, Definition 2.11, page 120]): 

The integer d is called the weight degree of / and denoted by wdeg(/). 

An ideal / C 1C[?7,X] is called quasi-homogeneous (of type (D, 1)) if there exists a set of quasi- 
homogeneous generators. In this case, we denote by ]K[?7,X]rf the EC-vector space of quasi-homoge- 
neous polynomials of weight degree d, and denote the set K[U,X]dnl. 

Proposition 1. Let I C K[?7,X] be an ideal. Then the following statements are equivalent: 

1. there exists a set of quasi-homogeneous generators of I; 

2. the sets 1^ are subspaces o/lC[?7,X]j, and I = ® defied- 
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Proof. Let f\,...,fp be a set of quasi -homogeneous generators of / and g = L'^ff £ ^- We 
prove first tliat / = ©(/gpj/^/. Let g(^) be its quasi-homogeneous component of weight degree 
d {i.e. g^^' = 2^ Therefore g = g^^'. Since g G /, there exists /ii /i^ G 

wdeg,{t)=d 

K[U,X] such that g = E^i^i/i- Since products of quasi-homogeneous polynomials are also 
quasi-homogeneous, 

gW = ^/j(<'-«"leg(/,))^. 

belongs to Id (where denotes the quasi-homogeneous component of h, of weight 

degree d — wdeg(^)) and hence / is equal to 0jgN/d. 

Conversely, let / be an ideal such that / = ^d^f^Id- Thus, for allJ G N and for all g El, g^^^ 
belongs to the ideal /. Let /i, . . . ,/p be a set of generators of the ideal /. Then the finite set 

of quasi-homogeneous polynomials F = {f-^^ : I < i < p,d eN} also generates /. Hence, / is 
generated by quasi-homogeneous polynomials. 

□ 

If / is a quasi-homogeneous ideal, then its weighted Hilbert series wHS/(/) G is defined as 
follows: 

wHS7(0 = L ^imiMU,X]d/Idy. 

den 

3 Transferring properties from to J^r 

In this section, we prove that generic structuraljjroperties (such as the dimension, the structure of 
the leading monomial ideal,. . .) of the ideal J'r are the same as properties of the ideal &r where 
several generic forms have been added. Hence several classical properties of the determinantal ideal 
&r transfer to the ideal J^^- For instance, this technique permits to obtain expUcit forms of the Hilbert 
series of the ideal J^^. 

In the following, we denote by b and c the following sets of parameters: 

b = {bf^ I f GMon(D,)t),l <^<?im}; 

c = {cf]\\<i<n,\<j <m,\<l<nm}. 

Also, gi,...,g„m € K{b,c)[U,X] are generic quasi-homogeneous forms of type {D,l) and of 
weight degree D: 

teMon{D,k) l<i<n 
\<j<m 

We let S>r denote the ideal + (gi, • • • ,gnm) C K(b, c)[?7,X]. Here and subsequently, for a = 
{uij) G K ^ ° we denote by 9a the following evaluation morphism: 

<Pa: IK[o] — > E 

/(ai,i,...,a„,m) I — > • • • ,««,m) 

Also, for (b,c) G K"*"^^ " ''^"'"^ we denote by v/b,c the evaluation morphism: 

VAb,c: K[b,c] K 
/(b,c) ^ /(b,c) 
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By abuse of notation, we let (Pa{^r) (resp. '/'b,c(^r)) denote the ideal {uij — (Paifij)) C 

We call property a map from the set of ideals of K[?7,X] to {true, false}: 
^: ldeals(E[f/,X]) {true, false} . 

Definition 1. Let hea property. We say that is 

• J^r-generic if there exists a non-empty Zariski open subset O C 



such that 



(9a(^r)) 



true; 



■nm{{°-^+')+nm) 



such that 



• ^r-generic if there exists a non-empty Zariski open subset O C 

(b,c) G O ^ ^ (V^b,c i^r)) = true. 

The following lemma is the main result of this section: 
Lemma 1. A property is J'r-generic if and only if it is Sir-generic. 

Proof. To obtain a representation of 9a (^'•) ^ generic a as a specialization of J^r (and conversely), 
it is sufficient to perform a hnear combination of the generators. The point of this proof is to show 
that genericity is preserved during this hnear transform. 

In the sequel we denote by 2t, *8 and £ the following matrices (of respective sizes nm x (^ ^j^^) , 
nm X and nm x nm): 



X2 X^ 



21 = 



05 



(nm) (nm) 
^(1) ^(1) 



(nm) Anm) 



(nm) 



(nm) 



' b 



Therefore, we have 



\^n,m fn,m / 
\gnm I 



(nm) 
^n,m j 



I Ml,l\ 



4 ) 



= ld„ 




-21- 



/ 4 \ 



+ 55- 



X^-'X2 

V 4 J 
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In this proof, for a € IC""'^"^^*) (resp. b G K'™(''"«^*),c G K"'™'), the notation A (resp. B,C) 
stands for the evaluation of the matrix 21 (resp. 55, C) at a (resp. b,c). Also, we implicitly identify A 
with a (resp. B with b, C with c, 21 with o, *B with b, C with c). 

• Let ^ be a. ^^-generic property. Thus there exists a non-zero polynomial /ii(2l) G ]K[a] such 
thatif /ji(A) /Othen ^ f<Pa(^r)) = true. 



Let adj(£) denote the adjugate of C (i.e. adj(£) = det(C) • C Mn IK(c)). Consider the polyno- 
mial hi defined by hi{^X) = /zi(-adj((j:) • 03) G K[b,c]. The inequality det{€)hi{^X) / 

defines a non-empty Zariski open subset O CK ° ^ Let (B, C) G C? be an element in 
this set, then C is invertible since det(C) ^ 0. Let A be the matrix A = — adj(C) • B. Therefore 
the generators of the ideal q)^ ("^'") invertible linear combination of the generators of 

V^b,c i^r]- Consequently, (piiJ^r) = Vhx i^r]- Moreover, /21(A) = /2i(B,C) / implies that 



the polynomial hi is not identically 0. Therefore, 

V(b,c)GO,^(vAb,c (¥r)) = ^ (^(Pi(Jr)) = true, 
and hence =^ is a i^^-generic property. 

Conversely, consider a ^^-generic property Thus, there exists a non-zero polynomial 
/J2(*B,C) GK[b,c] such that if /j2(b,c) /Othen ^ (^Vb,c{^r)^ = true. Since ^ is ^^-generic, 
there exists (b,c) such that /j2(b,c)det(c) 7^0. Let be the polynomial /j2(b) =/z2(— C-5S,C). 
Since det(C) 7^ 0, the matrix C is invertible and /i2(— C^^ • B) = /j2(B,C) 7^ and hence the 
polynomial h2 is not identically 0. Moreover, if a G » ) is such that /12(A) 7^ 0, then 

/i2(-C- A,C) 7^ and thus ^ ( '/''-C A.c(^r) ) = true. Finally, Y-C-A,c{^r) = 9a(=A) since 



the generators of Y-C A,c{^r) are an invertible linear combination of that of (Pa{<^r) (the linear 
transformation being given by the invertible matrix C) and hence they generate the same ideal. 
Therefore, the property ^ is J^^-generic. 

□ 

In the sequel, ~< is an admissible monomial ordering (see e.g f7\ Chapter 2, §2, Definition 1]), 
and for any polynomial /, LM(/) denotes its leading monomial with respect to ^. If / is an ideal 
of K[U,X], K{a)[U,X], or K{b,c)[U,X], we let LM(/) denote the ideal generated by the leading 
monomials of the polynomials. 

By slight abuse of notation, if h and /2 are ideals ofK[U,X],K{a)[U,X],orK{b,c)[U,X] {h and 
I2 are not necessarily ideals of the same ring), we write LM(/i) = LM(/2) if the sets {LM(/) | / G /i} 
and {LM(/) | / G /2} are equal. 

Lemma 2. Let 3^-^ and ff*-^ he the properties defined by 

( true if LM{I) = IM (Jr) ; 
I false otherwise. 



true if\-U{I) = LU[% 
false otherwise. 



8 



Then (resp. ^^^) is a J^r-g^neric (resp. generic) property. 

Proof. We prove here that is J^r-generic (the proof for is similar). 

The outhne of this proof is the following: during the computation of a Grobner basis G of 
in lC(a)[?7,X] (for instance with Buchberger's algorithm), a finite number of polynomials are con- 
structed. Let <Pa be a specialization. If the images by (pa of the leading coefficients of all non-zero 
polynomials arising during the computation do not vanish, then (Pa{G) C <Pa(=^r) is a Grobner basis 
of the ideal it generates. It remains to prove that (Pa{G) is a Grobner basis of <Pa(=^r)- This is achieved 
by showing that generically, the normal form (with respect to (Pa{G)) of the generators of (Pa{J^r) is 
equal to zero. 

For polynomials /i ,/2, we let LC(/i ) (resp. LC(/2)) denote the leading coefficient of fi (resp. /2) 

jc 1/^ ^\ LCM(LM(/i).LM(ft)) r LCM(LM(/i ),LM(/,)) ^ . <-u c i ■ i a r 

and Spol(/i ,/2) = — lc(/i) lm(/i) /i mMW)f^ A&notQ the S-polynomial of /i and /2. 

We need to prove that there exists a non-empty Zariski open subset Oi C IK ^ ' such that 

a G Oi ^ LM(9a(?r)) = LM(^,). 

To do so, consider a Grobner basis G C ]K(a) [?7,X] of J^,- such that each polynomial g can be written 
as a combination g = Y^hifi, where the //'s range over the set of minors of size r + 1 of ^ and 
the polynomials Uij — fij, and hg G ]fC[a][?7,X]. Buchberger's criterion states that S-polynomials of 
polynomials in a Grobner basis reduce to zero ||7] Chapter 2, §6, Theorem 6]. Thus each S-polynomial 
of gi,gj G G can be rewritten as an algebraic combination 

Spo\{gi,gj) = Y,f^'egi, 

i 

where the polynomials h'g belongs to ]K(o) [U,X] and such that i , . . . , g,, ^ } C G and for each 1 < 5 < 
tij, lM{gs) divides lM{Spo\{g,g') -Yj'el\h'ggi)- Next, consider: 

• the product Qi (a) = YlgeG '-C(g) of the leading coefficients of the polynomials in the Grobner 
basis; 

• for all {gi,gj) G such that Spo\{gi,gj) / 0, the product Q2{a) of the numerators and denom- 
inators of the leading coefficients arising during the reduction of Spo\{gi,g j). 

These coefficients belongs to lC[a]. Denote by 2(a) = Qi{a)Q2{a) G K[a\ their product. The 
inequality / defines a non-empty Zariski open subset Oi C K ^ ° ^IfaGOi,then 

t 

<Pa(Spol(g,/)) = ^ (pa{h'(.)(paigi), 

l=\ 

and for each \<i<t, LM((pa(g,)) divides LM(<Pa(Spolfe,/))^i:^-\ 9a(/j£)<Pafe))- Thus <Pa(G) is 
a Grobner basis of the ideal it spans. Moreover, {(pa{G)) C 9a (-^r)- 

We prove now that there exists a non-empty Zariski open set where the other inclusion (pa{^r) C 
{(Pa{G)) holds. Let NFg(-) be the normal form associated to this Grobner basis (as defined as the 
remainder of the division by G Chapter 2, §6, Proposition 1]). For each generator / of 

(i.e. either a maximal minor of the matrix or a polynomial uij — fj), we have that NFg(/) = 0. 
During the computation of N Fg(/) by using the division Algorithm in Q Chapter 2, §3], a finite set of 
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polynomials (in ]K(a)[?7,X]) is constructed. Let Qi, £ K[a\ denote the product of the numerators and 
denominators of all their nonzero coefficients. Consequently, if {a) ^ 0, then ^^(p^{G)i9a{f)) = 

and hence %(/) S {(pa{G)). Repeating this operation for all the generators of yields a finite set 

( f) 

of non-identically null polynomials GK[a]. Let Q4 G K[a] denote their product. Therefore, if 
e4(a) / 0, then (p^{Jr) C ((Pa(G)). 

Finally, consider the non-empty Zariski open subset O C 1C"'"( » ) defined by the inequality 
gi ■ 62 • 24 7^ 0. For all a G (9, we have (pa(=^r) = (<Pa(G))- 

□ 

Corollary 1. The leading monomials of J^r cif^ the same as that of 

LM f = LM 



Proof. By Lemmas 1 and 2 the property (resp. is ^^-generic and i^^-generic. Since 

(resp. is ^,--generic, there exists a non-empty Zariski open subset C?i C K ' ' 

(resp. O2 C such that, for (b,c) G Oi (resp. O2), LM (i/^(b,c)(^)) = LM 



(resp. LM = LM 

Notice that C?i n O2 is not empty, since for the Zariski topology, the intersection of finitely-many 
non-empty open subsets is non-empty. Let (b,c) be an element of Oi n C?2- Then 

LM (Jr) = LM (vA(b,c)(^-)) = LM (¥r) . 

□ 

Corollary 2. The weighted Hilbert series of ■J'r is the same as that of Sir- 
Proof. It is known that, for any positively graded ideal / and for any monomial ordering, wHS/(f) = 
wHSlm(/)(0 (s^^ ^-g- the proof of 17^ Chapter 9, §3, Proposition 9] which is also valid for quasi- 
homogeneous ideals). By Corollary jlj LM {^^r^ = LM {^^r^> which implies that 

wHS, ^,(~\ it) = wHS, (0, 
and hence wHS-^^ (t) = wHS^ (f). □ 



4 The case ^ > {n — r){m — r) 

As we will see in the sequel, the KruU dimension of the ring ]K(o)[X]/^^ is equal to max(A: — {n — 
r){m — r),0). This section is devoted to the study of the case k> {n — r){m — r). 

We show here that the algebraic structure of the ideal J^^ is closely related to that of a generic 
section of a determinantal variety. 

We recall that the polynomials gi are defined by 

si= I bfw £ eg.,,. 

teMon(D,k) l<i<n 
l<i<m 
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Lemma 3. Let \ <l< nm be an integer. If the polynomial gg divides zero in the ring 

K{b,c)[U,X]/{^^r + {gu...,gi-i)), 
then there exists a prime ideal P associated to {gi, - ■ ■ ,ge-i) such that dim(P) = 0. 

Proof. If g( divides zero in K{b,c)[U ,X]/ - + (^i, • • • then there exists a prime ideal P 

associated to + {gi, - ■ ■ ,ge-i) such that g( G P. For £ < nm, let b^-^''^ and c^-^-* denote the sets of 
parameters 

b(<^) = {bl'^ \t £Mon{D,k),\<s<£} 

c(^^) = {c\'^j \ I <i<n,l < j <m,\ <s<i}. 

Since + {gi ■ ^ge-i)) is an ideal of K{b^^'^-'^\c^^^-'^^)[U,X], and P is an associated prime, 
there exists a Grobner basis Gp of P (for any monomial ordering ^) which is a finite subset of 
K(b(^^-i),c(^^-i))[[/,X]. 

Let NFp(-) denote the normal form associated to this Grobner basis (as defined as the remainder 
of the division by Gp in Q Chapter 2, §6, Proposition 1]). 

Since g^ G P, we have H^p{ge) = 0. By linearity of NF/3(-), we obtain 

bf'NFp(0+ i: cS5nFp(m,j)=0. 

teMoa(D,k) l<i<n 
l<j<m 

Since Gp C K{b^-^^^\c^-^^^^)[U ,X], we can deduce that for any monomial t, 

Ufp{t)eK{b'^-^-^\c^^^'-^'>)[U,X]. 

Therefore, by algebraic independence of the parameters, the following properties hold: for all t G 
Mon(D,fc), NFp(f) = 0, and for all i,j, NFp(m,j) = 0. Consequently, all monomials of weight degree 
D in K(b, c) [U,X] are in P, and hence P has dimension 0. □ 

Lemma 4. For all i G {2, . . . , nm}, the polynomial gg does not divide zero in the ring 

K{b,c)[U,X]/i&r + {gi,-.-,ge-i)) 
and dim{% + {gi,--- ,g('}) =k + {n + m — r)r — £. 

Proof. We prove the Lemma by induction on £. According to |[T9l Corollary 2 of Theorem 1], 
the ring K{b,c)[U,X]/&,. is Cohen-Macaulay and purely equidimensional. First, notice that the 
dimension is equal to k + {n + m — r)r foT £ = since the dimension of the ideal C K[U] is 
{n + m — r)r (see e.g. fSl and references therein). Now, suppose that the dimension of the ideal 
% + {g\,---.gi>-\) C ]K(b,c)[?7,X] IS k + {n + m - r)r - £ + \. Since the ring ^{b,c)[U ,X]/ % is 
Cohen-Macaulay and (gi , . . . ,gt-\) has co-dimension £—\ in lC(b, c) [i/jX], the Macaulay unmixed- 
ness Theorem ||8l Corollary 18.14] implies that (gi, . . . has no embedded component and is 

equidimensional in ]K(b, c)[?7,X]/^,-. Hence + • • • as an ideal in ]K(b, c)[?7,X] has 

no embedded component and is equidimensional. By contradiction, suppose that g^ divides zero 
in ]K(b, c)[[/,X]/(^r + (gi, • • • By Lemma |3j there exists a prime P associated to + 

(gi, . . . such that dim(f ) = 0, which contradicts the fact that + (gi, . . . ,gi-\) is purely equidi- 

mensional of dimension k + {n + m — r)r — £ + \ > 0. □ 
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Lemma 5. The Hilbert series of the J'r C ]K(a) [X] equals the weighted Hilbert series of J'r C 
K(o)[Z,f/]. 

Proof. Let ^/gjf denote a lexicographical ordering on ]K(a)[X,?7] such that a;^. ^/g^t- m, j for all By 
E Section 9.3, Proposition 9], HS^r (f) = HSlm^, and ^WS j{t) = wHS|^|y^ Since 

LM^,^^(M,- y — fi,j) = Uij, we deduce that all monomials which are multiples of a variable y are in 
LM^i^^{^i.)- Therefore, the remaining monomials in LM^,^^.(=/r) are in ]K(a)[X]: 



lM^,Mr) = [{uij}UlM^,J,J^rnK{a)[X]) 
= <K,}ULM^,„(^,)). 



Therefore, isomorphic (as a graded ]K( a) -algebra) to l]^^^^^^- Thus 



and hence 

HS^^(0 = wHSj^(0. 



□ 



In the sequel, Ar{t) denotes the r x r matrix whose (/,7)-entry is Y,k {"\ ') {" t'')^'^- following 
theorem is the main result of this section: 

Theorem 1. The dimension of the ideal J^y is k — (n — r){m — r) and its Hilbert series is 

det(A,(f^))(l-f^)(«-'")('"-'-) 



t^iO{\-tf 



Proof. According to [6, Corollary 1] (and references therein), the ideal seen as an ideal of K[U] 
has dimension {m + n — r)r and its Hilbert series (for the standard gradation: deg(M, y) = 1) is the 
power series expansion of 

By putting a weight D on each variable Uij (i.e. deg(M,,y) = D), the weighted Hilbert series of &r C 
K[U] is 

uc r A detA,(f^) 

By considering as an ideal of c)[?7,X], the dimension becomes k + {m + n — r)r and its 
weighted Hilbert series is 

detA,-(f^) 

wH5^ c:K(b,c)[j/,x](0 - , : r- 

^ ^ tH2)(\-tf(\-tDyn+m-r)r 

According to Lemma |4j for each £ < nm, the polynomial gf does not divide zero in the ring 
K{a)[U,X]/{^^r + {gi,--- ,ge-i))- This implies the following relations: 

dim{% + {gu...,ge)) = dim(^, + (gi,...,g^_i))-l 
wHS5^,.+(g,,...,g^,)(0 = (l-?^)wHS5^,.+(^,,...^g^_,)(0. 
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Therefore the dimension of is ^ — nm +{n + m — r)r and its quasi-homogeneous Hilbert series is 

det(A,-(?^)) det(A,(?^))(l-?^)(«-'-)('"-'-) 

Wnb^U) = — r = -r-r . 

By Corollary [ijjhe ideal J^r has the same weightedHilbert series. Finally, by Lemma [s] the Hilbert 
series of J^r = ^r^) K{a) [X] is the same as that of J^r- □ 

Corollary 3. The degree of the ideal J'y is: 

DEG(^,) = d(«-'-) Y[ 



y^f^ {m-\-i)\{n-r + i)\ 



m-r-l (n+m-r-l\ 
_ jy{n-r)(m-r) T~i V r+i I 

'•=0 I j 

Proof. From |[T6l Example 14.4.14], the degree of the ideal is 



.^^ (m - 1 - /) ! (n - r + /) ! ' 



Since the degree is equal to the numerator of the Hilbert series of Qi,- evaluated at f = 1, 

detA,(l)= II 



;=0 {m - I - i)\{n - r + i)\' 
By Theorem [TJ the Hilbert series of J^^ is 

det(A,(?^)) 



detA,(f^)(l + f + • • . + fO-l)nm-(n+m-r)r 



t^iv ^1 f'^k—{n+m—r)r+nm 

Thus, the evaluation of the numerator in f = 1 yields 

deg(^,)=d("-'-)('"-'-) n ^ ^ 

i=0 



m — 1 — /) ! (n — r + /) ! 



To prove the second equality, notice that 

i\{n + m — r — i — 1) 

i=0 I ) i=o 

By substituting / by m — r — 1 — /, we obtain that 



fr ^ ^ '"n ' /!(n + m-r-/-l)! 

11 ^«+m-r-h 11 Cr + /)!(n + m-2r-/- 1)!' 



m—r—l m—r—l 

n {n + m-r- 

m—r—l m—r—l 

n = n i^-i-^v- 

m—r—l m— ; — 1 

Y\ {n + m-lr-i-\)\ = Y\ {n-r + i)\. 
i=0 !=0 
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Consequently, 



t=o {m-l-iy.{n-r + i)\ ^n+m-r-iy 

□ 



5 The over-determined case 

To study the over-determined case (k < {n — r){m — r)), we need to assume a variant of Froberg's 
conjecture flSl : 

Conjecture 1. Let denote the vector space of quasi-homogeneous polynomials of weight degree 
i in {gi,- ■ - jge)- Then the linear map 

K(b,c)[[/,X];/% K{b,t)[U,X]i+D/&t,i+D 

f ' — > fgi+\ 
has maximal rank, i.e. it is either injective or onto. 

Remark 1. Ifk + {n + m — r)r — ^ > 0, then Conjecture is proved by Lemma^ gf^i does not divide 
zero in ]K(b, c) [U,X]/ [Sir + {g\,... ,ge)) and hence the linear map is injective for all i G N. 

Notation. Given a power series S{t) € we let denote the power series obtained by 

truncated S{t) at its first non positive coefficient. 

Lemma 6. If Conjecture [7] is true, then the Hilbert series of + {gi^ - ■ ■ ,gi+\) is 

Proof. In this proof, for simplicity of notation, we let R denote the ring ]K(b,c)[?7,X]. If S{t) = 
Y^ien^i^' £ ^[M] is a power series, [5'(f)]>Q denotes the series 

[5(0]>o = Imax(.,-,Oy. 

Let ann(g^+i) be the ideal {f eR: fge+i & ^r+{gi,---,ge)}- For / G N, consider the following exact 
sequence: 

^ ann(^^+i),- ^ /?,•/% Ri+o/^u+D ^ 

^Ri+o/^i+ij+D^O. 

Conjecture [T] states that 

dim(ann(g^+i)/) = max(0,dim(/?,-/^/,;) -dim(/?,■+D/^f,!■+D))• 
Since the alternate sum of the dimensions of the vector spaces occurring in an exact sequence is zero, 
it follows that 

dim{Ri+D/^i+\,i+D) = dim{Ri+D/&ej+D)-dim{Ri/&e,i)+ 

max(0,dim(/?,-/%,-) - dim(/?,-+z)/%,-+z))) 
= max(0,dim(/?,-+z)/^f,,-+z)) -dim(/?i/^^,-))- 
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Multiplying this equation by t'^ yields 

[?'+^]wHS^^+(g,,...,g,^,)(0 = dim{Ri+Dl&t+u+D)) 

= max (0, dim(/?;+z)/^f,,-+z)) - dim{Ri/&ej)) 
= max (0, [f'+^](l -f^) wHS^^+(,.,...,,,)(0) 

= r^][(i-^^)wHs^,^+^,„...,,,)(o]>o. 

Since any monomial in K{a)[X,U] of weight degree greater that D is a multiple of a monomial of 
weight degree D, we deduce that if there exists iq > D such that 

[f'»]wHS5,,+(,„...,,,,,)(0 = 0, 
then for all / > iq, [?'] wHS^^^+^gi ^^^^^ (?) = 0. Therefore 

[t'^^] wHS^,.+^,„...,,,,^.)(0 = [(1 -?^)wHS^^+fe„...,,,,)(0]^ , 
Finally, by summing over /, we get 



□ 



Theorem 2. TfConjecfMrejiJ/s frwe, the Hilbert series of is 

HS,y,(0 



where Ar{t) is the rxr matrix whose {i,j)-entry is ^ 



f^u)(l-0*^ 

min(m— i,n— 7) 



k=Q 



m-i\ (n-]\^^. 



Proof. By applying nm times Lemma |6] we obtain that 



wHS^^(0 



detA,(f°) 



Let 5 = ^o<! ^i^' S ^[M] be a power series such that qq > 0, and let /q G N U {00} be defined as 



00 if for all / > 0, a; > 0; 
min({/ 1 Ui < 0}) otherwise. 



Therefore, = Lo<!<!o'^f^'- convention, for / < 0, we put a, = 0. Then 

(l-f^)5(0 = ^o<iiai-ai-D)t' 

(i-f^)[5(0]+ = i:o<i<iMi-^i-Dy ■ 

Consequently, the coefficients of (1 — t^)S{t) and of (1 — t^) [S{t)]_^_ are equal up to the index /q- 
• If /o = 00, then (1 = (1 -f^) [5(0]+ and hence 

[(1-^^)5(0], =[(1- [5(0]+], ; 
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• if /() < °°, then ai^^o is positive and thus — ai^^o is negative. Let ii be the index of the first 
non-positive coefficient of (1 — t^)S{t). Then ii < iq, and hence 

[(i-f^)5(o]+ = [(i-^^)[5(ou+. 

Therefore, for all power series S G such that 5(0) > 0, we have 

[(1 -?^) [5] [(1-^^)5]^. 
Consequently, an induction shows that 



wHS^(0 



:i-t 



D\ (n—r) (m—r) 



r) detA(?«) 



Then, by Corollaryj2| wHS^ (?) = wHS-^(f). Finally, by Lemmajsj we conclude that HS,j?;(f) = 



6 Complexity analysis 

Using the previous results on the Hilbert series of .^r, we analyze now the arithmetic complexity of 
solving the generalized MinRank problem with Grobner bases algorithms. In the first part of this 



section (until Section 6.2 1, we consider the homogeneous MinRank problem (i.e. the polynomials 
are homogeneous). 

Computing a Grobner basis of the ideal (Pai^r) for the lexicographical ordering yields an explicit 
description of the set of points V such that the matrix 



/<Pa(/l,l) ••• <Pa(/l 



<Pa 



1 ,m I 



V<Pa(/n,l) ••• %(/«,, 



has rank less than r + 1 . In this section, we study the complexity of this computation when a € W" 

— 11 in 

is generic (i.e. a belongs to a given non-empty Zariski open subset of K ) by using the theoretical 
results from Sections |4] and [5] We focus on the 0-dimensional cases k = {n — r){m — r) and k < 
{n — r){m — r) (over-determined case). Therefore, the set of points where the evaluation of the matrix 
has rank less than r + 1 is finite. 
In order to compute this set of points, we use the following strategy: 

• compute a Grobner basis of <Pa(-^r) for the grevlex (graded reverse lexicographical) ordering 
with the 7^5 algorithm ||9l ; 

• convert it into a lexicographical Grobner basis of 9a(=^r) by using the FGLM algorithm lfT0l[T4]| . 

First, we recall some results about the complexity of the algorithms F5 and FGLM. The two 
quantities which allow us to estimate their complexity are respectively the degree of regularity and 
the degree of the ideal. The degree of regularity of a homogeneous ideal I is the maximum degree 
in a minimal Grobner basis (and hence in general it depends on the monomial ordering). If 7 is a 
0-dimensional ideal, it is also the smallest integer d such that all monomials of degree are in 7 and it 
is independent on the monomial ordering. Moreover, in the 0-dimensional case, the Hilbert series is a 
polynomial from which the degree of regularity can be read off: 



"regV 



:deg(HS/(0) + l. 
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In the sequel, ft) denotes a feasible exponent for the matrix multiplication (i.e. a number such that 
there exists an deterministic algorithm which computes the product of two n xn matrices in 0{n'^) 
arithmetic operations in IC). 

The following proposition and its proof are a variant of a result known in the context of semi- 
regular sequences (see e.g. [1 , Proposition 10]): 

Proposition 2 dUD- Let hi,. . .,h( £ K[xi ,Xk] be homogeneous polynomials of degrees di,. . .,d£, 
and I = {hi, . . . ,hf). The complexity of computing a Grobner basis of I for a monomial ordering -< is 
upper bounded by 

^ ( / A+Preg(/)\ nprrnV"' A+ID'reg(/)\ ^ A+Preg(/)-^A\ 

""[[y D.eg(/) (, D.g(/) )hX i^.g(/)-^,- ))■ 

Proof. Since / is homogeneous, a Grobner basis can be obtained by computing the row echelon form 
of the so-called Macaulay matrix of the system up to degree Dreg(^)- This matrix is constructed as 
follows: 

• the rows are indexed by the products thi, where I <i <t and t G 'K[xi , . . . ,;ci:] is a monomial of 
degree at most ID>reg(^) — df, 

• the columns are indexed by the monomials m € lC[xi, . . . ,Xi-] of degree at most ID)reg(^) and are 
sorted in descending order with respect to -< ; 

• the coefficient at the intersection of the row thi and the column m is the coefficient of m in the 
polynomial thj. 

The number of columns of this matrix is the number of monomials in ]IC[xi, . . . ,x^] of degree at 
most Dreg(^)> namely {^^'"^jI^)- The number of rows is ^f^j and its rank is equal to 

((X:^')-DEG(/)). "" 

According to ll24l Theorem 2.10], the complexity of computing the row echelon form of a. p xq 
matrix of rank r is upper bounded by 0{r'^^^pq). 

Consequently, the complexity of computing a Grobner basis of / is upper bounded by 

""[[y D.eg(/) i D,eg(/) )hX ^r.,{I)-di j)' 

□ 

Remark 2. Notice that 

V ID),eg(/) J V ID)reg(/) J 

' h + D,,^iI)-dA ^ A + D,eg(/)\ 

ti\ o,,,ii)-dt J - \ Dreg(/) 

Therefore, the complexity of computing a Grobner basis of I can also be upper bounded by the simpler 



expression O 
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Lemma 7. If k = [n — r){m — r), then the degree of regularity of J'r is 

Dreg Dr(m - r) + (D - 1 )/t + 1 . 

Proof. According to Theorem [TJ the Hilbert series of J^r is 

detA,(f^)(l-f^)("-'')('«-'") 



HS,y,(0 



By definition of the matrix the highest degree on each row is reached on the diagonal. Thus, 

the degree of det(A,.(f)) is the degree of the product of its diagonal elements: 



deg(det(Ar(f))) = ^(min(?i,wi) — /) = rrn 

!=1 



r+1 
2 



Therefore, we can compute the degree of the Hilbert series which is a polynomial since the ideal is 
0-dimensional: 

Oreg(=A) = deg(HS^/^(0) + l 

= deg(det(A,(f^)))+D(?i-r)(m-r)-D0 -k+l 
= D(rm- +nm-{n + m-r)r- {^^)-k+l 
= Dr{m-r) + {D-\)k+\. 

□ 

Corollary 4. lfk={n — r){m — r), then there exists a non-empty Zariski open subset C? C K ^ ' 
such that for all a G O, the degree of regularity of(pa{J^r) is 

Oreg {(p^{y,-))=Dr{m-r) + {D-l)k + \. 

Proof. According to Lemma|2j there exists a Zariski open subset O such that for all a G O, LM ( J^^) = 
LM {(pa{^r))- Consequently, the polynomials in minimal Grobner bases of and (Pa{^r) have the 
same leading monomials. Since the degree of regularity is the highest degree of the polynomials in a 
minimal Grobner basis, we have Dreg {(pai^r)) = ID'reg i^r)- Lemma |7] concludes the proof. □ 

The degree of regularity governs the complexity of the Grobner basis computation with respect 
to the grevlex ordering. The complexity of the algorithm FGLM is upper bounded by 0{k • DEG(/)^) 
which is polynomial in the degree of the ideal |[T0llT4ll . 

We can now state the main complexity result: 

Theorem 3. There exists a non-empty Zariski open subset O C 1C"'"( " ) such that for any a G O, 
the arithmetic complexity of computing a lexicographical Grobner basis of the ideal generated by the 
(r + 1) X (r + \)-minors of the matrix (pa{^) is upper bounded by 

where 2 < (O <3 is a feasible exponent for the matrix multiplication, and 
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ifk= {n — r){m — r), then 

Dreg = deg(HS^^(^y,.) (0) + 1 = Dr{m -r) + {D-\)k + l 

and DEG(<Pa(^.)) = HS,,(^,^)(1) = D— ("+'"-)^n;ro^-^ (nr-SS^r+i)r 
ifk < {n — r){m — r), then assuming that Conjecture^is true, 

D,eg = deg(HS^(^^)(0) + l 
and DEG(<Pa(-^r)) = HS^^(./^)(1) where 



D\nm-(n+m-r)r detA(f^) 



fO(0(l_f)<: 



Proof. The number of (r+ l)-minors of the matrix (pa(^) is Consequently, the theorem 

is a straightforward consequence of the bounds on the complexity of the F5 algorithm (Proposition [2]l 
and of the FGLM algorithm ifTOlfT?! . together with the formulas for the degree of regularity (Corollary 
[4]) and for the degree (Corollary [3]l. □ 

Remark 3. We can find a polynomial h{a) in Z[a] when the characteristic o/K is 0, such that 

/i(a) / O^aG O. 

Also note that this polynomial does not depend on the field K; ;/ IC = is a finite field (q = j, 
then the polynomial h{a) (where all coefficients are taken modulo p) verifies the requested property. 

Schwartz-Zippel's Lemma states that, if a is chosen uniformly at random in , the probability 

that h{a) =0 is upper bounded by deg{h)/q and therefore tends towards when the cardinality q 
of the field tends to infinity. This explains why these complexity results can be used for practical 
applications when char(]K) =0 orKis a sufficiently large finite field. 

6.1 Positive dimension 

When k > {n — r){m — r), the ideal .J^r has positive dimension. To achieve complexity bounds in that 
case, we need upper bounds on the degree of regularity . 

Lemma 8. If k > [n — r){m — r), then the degree of regularity of is 

Dreg(=/r) = Dr(m - r) + (D - 1 ) (« - r) (m - r) + 1 . 

Proof. Consider the ideal J obtained by specializing the last k— {n — r){m — r) to zero in J^^- We prove 
now that LM( J^^) = LM(7). First, notice that for the grevlex ordering, LM(7) C LM(J^,-). According 
to Theoremjlj the Hilbert series of the ideal J n K{a) [xi,... ,X(„-r)(m-?-)] equal to 

detA,(?^)(l-f^)("-'-)('"-'') 
^0(2) ^2 -?)("-'•)("'-'") 

By construction, J C K{a) [xi , . . . ,X(„_^)(„_^)], thus the Hilbert series of J as an ideal of ]K(o) [xi , . . . ,xt] 
is equal to 

detA,(f^)(l-f°)("-'-)("'-'-) 
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which is equal to the Hilbert series of ,^r- 

Since WS](t) = HSy,.(0 and LM(7) C LM(J^,-), we can deduce that LM(7) = LM(j?;). Con- 
sequently, the leading monomials in minimal Grobner bases of J and J^,- are the same. Hence, the 
polynomials in both Grobner bases have the same degrees since they are homogeneous. 

Finally, notice that the Grobner basis of the ideal J is the same as that of 7n]K(a) [xi , . . . ,X(„_r)(,„_^)] 
which, by Lemma|7| is a zero-dimensional ideal whose degree of regularity is Dr(m — r) + (D — 1 )(« — 
r) (m — r) + 1 . Therefore the maximal degree of the polynomials in the minimal reduced Grobner basis 
of J^r is Dr{m - r) + {D - \){n - r){m - r) + \. □ 

Using exactly the same argumentation as in the proof of Corollary [4| we deduce that 



Corollary 5. Ifk>{n — r){m — r), then there exists a non-empty Zariski open subset O C 
such that, for a G O, the maximal degree of the polynomials in a minimal grevlex Grobner basis of 

(Pa{Jr) is 

B)reg (9a(=^r)) = Dr{m - r) + {D - l){n - r){m - r) + I. 

Theorem 4. Ifk> {n — r){m — r), then there exists a non-empty Zariski open subset O C K""'^ ^ ) 
such that for any a & O, the arithmetic complexity of computing a grevlex Grobner basis of (p^i'^r) is 
upper bounded by 

/ / n \ / m \ fDr{m-r) + {D-\){n-r){m-r) + l+k 
^[\r+\)\r+\)\ k 

Proof. This is a consequence of Proposition [2] and Corollary [5] □ 
6.2 The 0-dimensional affine case 

For practical applications, the affine case (i.e. when the entries of the input matrix ^ are affine 
polynomials of degree D) is more often encountered than the homogeneous one. In this case, the 
matrix ^ is defined as follows 



l,m 



D 



t. 



\/n,l ••• fn 



^=OfeMon(f,/t) 



We show in this section that the complexity results (Theorems [3] and [4]l still hold in the affine case. 
This is achieved by considering the homogenized system: 

Definition 2. [7. Chapter 8, §2, Proposition 7] Let {qi,..., q^) £ lC[xi , . . . jX/^Y be an affine polynomial 
system. We let {qi, . . . ,qi/) G IK[xi , . . . ,Xk,Xk^iY denote its homogenized system defined by 

V/, s.t. l<i< £,qi{xi,...,Xk,Xk+i) =xff^'^'^qi ( ] . 

Notice that if an affine polynomial system has solutions, then the dimension of the ideal generated 
by its homogenized system is positive. 

The study of the homogenized system is motivated by the fact that, for the grevlex ordering, the 
dehomogenization of a Grobner basis of ((71, ... is a Grobner basis of {qi,... ,qe). Therefore, in 
order to compute a Grobner basis of the affine system, it is sufficient to compute a Grobner basis of 
the homogenized system (for which we have complexity estimates by Theorems |3] and |4]l. 

To estimate the complexity of the change of ordering, we need bounds on the degree of the ideal 
in the affine case: 
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Lemma 9. The degree of the ideal {qi,..., q^) is upper bounded by that of{qi,..., qi). 

Proof. The rings K[xi,. . . ,Xk]/{qu. ■ .,qi) and K[xi,. . . ,Xk,Xk+i]/{qi,. ■ . ,qe,Xk+i - 1) are isomor- 
phic. Therefore the degrees of the ideals {qi,...,q(i) and {qi,...,qg,Xk+i — 1) are equal. Since 
deg{xic+i — 1) = 1, we obtain: 

DEG((^i,...,^^)) = DEG{{qi,...,qi,Xk+i-l)) 
< DEGi{qu...,qi))- 

□ 

Lemma 10. The degree of regularity with respect to the grevlex ordering of the ideal {qi,-.. ,qi) is 
upper bounded by that of {q\,. . ., qi). 

Proof Let X denote the dehomogenization morphism: 

X: K[xi, . . . ,Xk+i] — > K[xi,...,Xk] 

f{xi,...,Xk,Xk+l) I — > f{xi,...,Xk,l) 

If G is a grevlex Grobner basis of (<^i , . . . , ^) , then x (G) is a grevlex Grobner basis of (^i , . . . , q() (this 
is a consequence of the following property of the grevlex ordering: V/ G ]K[;ici, . . . ,Xk+i] homogeneous, 
^^ixif)) = /C('-M(/))). Also, notice that for each g EG, any relation g = J^-^j qihi gives a relation 
Xis) =l4=iX{<li)Xihi) of lower degree since 

deg{x{qi)x{hi)) <deg{qihi). 

Consequently, a Grobner basis of (^i, . . . ,^^) can be obtained by computing the row echelon form of 
the Macaulay matrix of (^i, . . . ,qe) in degree Dreg((^i) • Therefore, the degree of regularity 

with respect to the grevlex ordering of the ideal (^i, . . . ,qi) is upper bounded by that of (^i, . . . ,^). 

□ 

We can now state the main complexity result for the affine generalized MinRank problem: 
Theorem 5. Suppose that the matrix ^ contains generic affine polynomials of degree D: 



\fn,l ■ ■ ■ fn,m/ 



=OteMoa(e,k) 



There exists a non identically null polynomial h G ]fC[o] such that for any a G » ) such that 
h{a) 7^ 0, the overall arithmetic complexity of computing the set of points such that the matrix (PaC-^) 
has rank less than r+l with Grobner basis algorithms is upper bounded by 

<C:,)C:,)(""rT-(-(-<-'))' 

where 2 < CO < 3 is a feasible exponent for the matrix multiplication and 
• ifk = {n — r){m — r), then 

Brsg<Dr{m-r) + {D-\)k+l, 

DEG(^(^,))<D("-')("-) n („_Jr;^:,^t,y 
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ifk < {n — r){m — r), then assuming that Conjecture^is true, 

lD)reg<deg(/'(0) + l, 

and DEG(<Pa(^r)) < ^(1) where 



Pit) 



(I _ fD\{n-r){m-r) detA(f^) 

^^(0(l-^)^ 



Proof. This is a direct consequence of Proposition [2j Lemma[9} Lemma 10 and the complexity of the 
FGLM algorithm HOlIll (0(A:DEG(^a(^r)))- □ 



7 Case studies 

The aim of this section is to compare the complexity of the grevlex Grobner basis computation with 
the degree of the ideal in the 0-dimensional case (i.e. the number of solutions of the MinRank problem 
counted with multiplicities). Since the "arithmetic" size (i.e. the number of coefficients) of the lexico- 
graphical Grobner basis is close to the degree of the ideal in the 0-dimensional case, it is interesting to 
identify families of parameters for which the arithmetic complexity of the computation is polynomial 
in this degree under genericity assumptions. 

Throughout this section, we focus on the 0-dimensional case: k= {n — r){m — r). Under genericity 
assumptions, we recall that, by Corollary [3] and Lemma|7} 



DEG 



Dr{m-r) + [D - \)k + \ 

f=o {m-\-i)\{n-r + i)V 



According to Theorem[5j the complexity of the computation of the grevlex Grobner basis is then upper 
bounded by 



O 



r+ 1 1 \r+ 1 



m \ (Dr{m-r) + {D-\)k+\ 



+ ^(DEG((pa(^.))y 



In this section, Q. and O are the Landau notations: for any positive functions / and g, we write 
f = Q.{g) (resp. / = 0{g)) if there exists a positive constant C such that f>C-g (resp. f <C-g). 



7.1 D grows, n, m, r are fixed 

We first study the case where n, m and r are fixed (and thus k = {n — r){m — r) is constant too), and D 
grows. In that case, the arithmetic complexity of the grevlex Grobner basis computation is 0(D*^®), 
and the degree is Q.{D'^). Therefore the arithmetic complexity has a polynomial dependence in the 
degree for these parameters. 



7.2 n grows, m, r, D are fixed 

This section is devoted to the study of the subfamilies of Generalized MinRank problems when the 
parameters m, r and D are constant values and n grows. Let i denote the constant value £ = m — r. 
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First, we assume that D = 1. When n grows, by Corollary [3] we have 

log(DEG) = log(^n|;^j 

r^log(«) 



n—^oo 



On the other hand. 



, ,^ . f (n-r)i + r£+l\ , f n \ , / m 

log(Compl) = cologf J+l°H''+lJ+ Hr+l 

, /«^+l\ , f n \ , / m 
= ^^°nr^+lj+^°nr+lj+^°nr+l 



(w(r£+l)+r+l)log(?i). 



Therefore, log(Compl) / log(DEG) ~ ^+ ) + — hence the number of arithmetic opera- 
tions is polynomial in the degree of the ideal. 

Also, if D > 2 is constant, a similar analysis yields 

log(DEG) = („-rKlog(D)+logfnj^^) 



, (=0 



17— >oo 



log(D)^n. 



. „ . fk + Drt + {D-\)k + l\ , f n \ , f m 

log(Compl) = colog( ^ ) +1°8L+1 ) +1°8L+1 

, /D?i^+1\ , ( n \ , / m 

= ^^°n(.-rKj+l°<.+ lJ+l°<.+ l 

, (Dni' 
~ CO log 

Then, using the fact that ( ) ~ « (alog(a) - j8 logfjS) - (a - j8) log(a - j8)), we obtain that 
log(Compl) ~ ?iw£(Dlog(D)-(D-l)log(D-l)). 

Therefore, log(Compl) / log(DEG) is upper bounded by a constant value and hence the arithmetic 
complexity of the Grobner basis computation is also polynomial in the degree of the ideal for this 
subclass of Generalized MinRank problems under genericity assumptions. 



7.3 The case r = m — 1 



The case r = m — 1 is a special case of the setting studied in Section 7.2 which arises in several 
applications, since it is the problem of finding at which points the evaluation of a polynomial matrix 
is rank defective. In this setting, the formulas in Theorem|5]are much simpler: 

• the 0-dimensional condition yields k = n — m + \; 

• ID^reg ^ — (n — m) ; 

• DEG<D"-'"+M " 
m — 1 
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(n,m,D,r,k) 


degree 


D,eg 


F4 time(Ma 


(6,5,2,4,2) 


60 


11 


0.001s 


(6,5,3,4,2) 


135 


17 


0.002s 


(6,5,4,4,2) 


240 


23 


0.004s 


(5,5,2,3,4) 


800 


17 


0.25s 


(8,5,2,4,4) 


1120 


13 


0.7s 


(5,5,3,3,4) 


4050 


27 


6.7s 


(6,5,2,3,6) 


11200 


19 


479s 



reg F4 time(Magma)FGLM time(Magma)F5 time/nb.ops(FGb)FGLM time(FGb) 



O.OOs/2'^-^" 
O.OOs/2 15-29 

0.01s/2i^-^'^ 
0.24S/225-56 

0.43S/226-71 

5.43s/23''-^^ 
94.85S/235.7 



0.001s 
0.019s 
0.09s 
6.3s 
20s 
567s 
17703s 



0.00s 
0.00s 
0.01s 
0.19s 
0.58s 

3s 
203s 



Table 1 : Experimental results 



Therefore, the arithmetic complexity of the Grobner basis computation is CompI = ). 
If D > 1 and m are fixed, log ( ( ^"^\)'^] ~ ft) log (^") and a direct application of Stirling's 

\^\n in-\-i/ J ^ — \ n / 

formula shows that 

cologf^") ~ a)(DlogD-(D-l)log(D-l))«. 

\ n J n^oo 

On the other hand, log(DEG) ~ nlogD. Therefore, log(Compl)/log(DEG) has a finite limit 
when n grows and m is fixed, showing that, in this setting, the arithmetic complexity is polynomial in 
the degree of the ideal. 

7.4 Experimental results 

In this section, we present some experimental results obtained by using the Grobner bases package 
FGb (using the F5 algorithm) and the implementation of the F4 algorithm in the MAGMA computer al- 
gebra system [3]. All instances were constructed as random (with uniform distribution) 0-dimensional 
MinRank problems (i.e. k = {n — r){m — r)) over the finite field F65521. All experiments were con- 
ducted on a 2.93 GHz Intel Xeon with 132 GB RAM. 

Useful information can be read from Table [T] First, the experimental values of the degree of 
regularity and of the degree match exactly the theoretical values given in Lemma [7] and in Corol- 
lary [3] Also, it can be noted that the most relevant indicator of the complexity of the Grobner basis 
computation seems to be the degree of the ideal. 

The comparison between the complexity bound and the degree of the ideal is illustrated in Fig- 
ures [T] and |2] First, Figure [T] shows that the bound on the complexity of the Grobner computa- 
tion is polynomial in the degree of the ideal when D grows (« = m = 20, r = 10 fixed), since 
log(ComplpJ/log(DEG) is upper bounded by 5. This is in accordance with the analysis performed 
in Section im 

Then Figure [2] shows empirically that if m = [j3?ij and r = [an\ — 1 (with a < j3 < 1) and n 
grows, then the complexity bound is also polynomial in the degree of the ideal. 

However, there also exist families of generalized MinRank problem where the complexity bound 
for the Grobner basis computation is not polynomial in the degree of ideal. For instance, taking n = m 
and fixing the values of r and D yields such a family. 

The experimental behavior of log(Complp^)/log(DEG) is plotted in Figure |3] We would like to 
point out that this does not necessarily mean that the complexity of the Grobner basis computation is 

(r+l) (r+l) ( k^^) ) 

not sharp and the figure only shows that the bound is not polynomial. 
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5 




bjO 



10 20 30 40 

D 

Figure 1: Numerical values of log(ComplpJ/log(DEG), for n = m = 20, r = I0,k = {n — r){m — r). 




Figure 2: Numerical values of log(ComplFj)/log(DEG), for m = [pn\,r = [an\ — l,D = l,k = 
(n — r){m — r). 
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Figure 3: Numerical values of log(Complpj)/log(DEG), for m = [j8«J ,r = \_an\ —l,D= l,k = 
(n — r){m — r). 



The problem of showing whether the actual arithmetic complexity of the F5 algorithm is poly- 
nomial or not in the degree of the ideal for any famiUes of parameters of the generalized MinRank 
problem remains an open problem. 

8 Application to bi-homogeneous systems of bi-degree (D, 1 ) 

In this section, we show that the previous complexity analysis can be used to obtain bounds on the 
complexity of solving bi-homogeneous systems of bi-degree (-0,1) by using Grobner bases algo- 
rithms. These structured systems can appear naturally in some applications, for instance in Geometry 
and in Optimization. Indeed the classical technique of Lagrange multipliers - when used to opti- 
mize a polynomial function under polynomial constraints - gives rise to a bi-homogeneous system of 
bi-degree {D,l). 

Bi-homogeneous polynomials are defined as follows: given two finite sets of variables X = 
{xo, ...,Xn^} and Y = {yo, ...,y„^}, a polynomial / G ]K[X,y] is called bi-homogeneous if for any 
A , G K, there exist dx,dy eN such that 

/(AX,mF) = AV^/(X,7). 

The couple {dx,dy) is called the bi-degree of /. 

In this section, we focus on generic systems of % -I- riy bi-homogeneous equations of bi-degree 
(Z), I). Such systems have a finite number of solutions on the biprojective space P"' x F'^y. One way 
to compute them is to start by computing their projection on P"% and then lift them to P""^ x P"^ 
by solving linear systems (this can be done since the equations are linear with respect the variables 

yo,---,yny)- 

The following proposition shows that computing the projection on F"y can be computed by solving 
a homogeneous MinRank problem. 
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Proposition 3. Let /i, . . . ,/m S IK[X,y] be a bi-homogeneous system of bi-degree (D, 1). Ifm > Uy, 
then (xo : ... : Xn^,yo : ... : y,^) E P"»' x P"> /i' a zero of this system if and only if the matrix 



jaCy (jtO) • • • -i^n^f 



dyo 


dfl \ 

dyny 


\ dyo ■ 


dfm 
dyny 1 



is rank defective. 
Proof. First, notice that 



\fm) 



:jacy(xo, 



\yny) 



Therefore, (xq : . . . : x„^.,3'o : • • • : 3'«J £ P"'' x P"' is a zero of the system if and only if (3^0, •. . ,3'n,,) 
belongs to the kernel of jacy . Since m> ny, the number of rows is greater than or equal to the number 
of columns of jacy , and hence jacy is rank defective. □ 

In applications, most of bi-homogeneous systems occurring are afftne: a polynomial / in the ring 
K[xi , . . . ,Xn^,yi is called affine of bi-degree (D, 1) if there exists a bi-homogeneous polyno- 
mial G K[xo, . . . ,Xn^.,yo, . . . , j„J of bi-degree (D, 1) such that 

f{xi,...,xn,,y\,...,yny) =f''{^,Xl,...,Xn,,l,yl,...,yny)■ 
This means that each monomial occurring in / has bi-degree {i,j) with / < D and 7 < 1. Notice that 
the polynomial /'^ is uniquely defined and that Proposition [s] also holds in the affine context: 

Proposition 4. Let /i,. ..,/,„ € K[xi, . . . ,x„^,y\, . . . ,yn^] be an affine system of bi-degree {D,l). If 
m > ny and {x[,. . . ,Xn^,yi ,.. . ,yn^,) £ M."-" x is a zero of the system, then the mx (ny + \) matrix 



jaCy [xi , . . . ,x„^ 



f fl{xi,...,Xn^,0,...,0) ^ 

\fm{xi,...,x„^,0,...,0) 



dU 

dyo 



dfl \ 

dy„y \ 



df„. 

dyny I 



is rank defective. 

Proof. The proof is similar to that of [3] since 



:jaCy(xi, 



\fmj 



/1\ 

y\ 

\yny) 



Therefore, if {x\,. . . ,Xn^,y\, . . . ,yn^) is a zero of the system then there is a non-zero vector in the 
kernel of jaCy (however in the affine case, the converse is not true). Since m > ny, the number of rows 
is greater than or equal to the number of columns of jaCy , and hence jaCy is rank defective. □ 
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An algebraic description of the variety V of a 0-dimensional polynomial system can be obtained 
by computing a rational parametrization, i.e. a polynomial g{u) G K[u\ and a set of rational functions 

gi,...,gn^,hi,...,hn^, €K{u) such that 

{xi,...,Xn^,yi,...,ynJ eV 

t 

3u £ K,s.t.g{u) =0,V/ e {l,...,nx},Xi = gi{u),yj £ {I, . . . ,ny},yj = hj{u). 

To obtain a rational parametrization, we need a separating element: a linear form which takes 
different values on all points of V. Therefore, a rational parametrization exists only if the cardinality 
of the field K is or large enough. 

Under the assumption that the field K is sufficiently large, Algorithm[T]uses the property described 
in Proposition [4] to find a rational parametrization of the zeroes of a radical and 0-dimensional system 
of rix + riy affine polynomials of bi-degree {D,l). The algorithm proceeds by first computing a rational 
parametrization of the projection of the zero set on K"' . This is done by computing a lexicographical 
Grobner basis of a Generalized MinRank Problem. Then this parametrization is lifted to the whole 
space by solving a linear system (this can be done since the equations are linear with respect to the 
variables J 1,..., 

The success of Algorithm [T] depends on the choice of the parameters a (a linear change of co- 
ordinates such that x„ is a separating element) and M. However, as we will see in Theorem |6] if the 
cardinality of K is infinite or large enough, then almost all choices of a and M are good. Therefore, 
these parameters can be chosen at random. If Algorithm [T] unluckily fails, then it can be restarted with 
the same algebraic system and different values of a and M. 
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Algorithm 1 Rational parametrization of systems of bi-degree (D, 1) 

Require: /i, . . . ^fn^+n^. S ]K[X,F] a system of affine polynomials of bi-degree (D, 1) such that the 
ideal they generate is radical and 0-dimensional; 
(ai,...,a„,._i) 

a full rank matrix M = {mtj) € lC">'>'("^+%). 
Ensure: Returns a rational parametrization of the variety of the system or "fail". 
1: Compute for each / G {1, . . . ,nx + ny}, 

fi{x\,...,Xn^-i,u,yi,...,yn,) = fi{xi,...,Xn^^i,u- ^ aexe,yi , . . . ,yny) ■ 

i=\ 

2: Compute the matrix jacy(7i , . . . 

3: Compute a lex Grobner basis G of the ideal / C "K^xi,. . . ,Xn^^\,u] generated by the maximal 

minors of the matrix jacy(/i, . . . ,/«_,+„ J. If the Grobner basis has the following shape (the shape 
position): 

xi -gi{u) 
xi-giiu) 

Xn^-\ - 8n,-\{u) 

then continue to Step 4, else return "fail". 
4: Using M, compute a linear combination of the polynomials of the system evaluated at 



/ fi{yi,---,yn,,u)' 



\fn,iyu---,yny,u)^ 



■ M- 



I fi{gi{u),---,gn,-i{u),u,yi,...,y„^.) modg{u) 
\/«,+«v(§i ("),••• ,gn,-i {u),u,yi ,yn,) mod g{u)y 



5: If the linear system f\ = ... = f„^_ = has rank Uy (as a linear system in K{u) [Y] where the 

variables are 3^1 , . . . continue to Step 6, else return "fail". 
6: Using Cramer's rule, solve the system /i = . . . = = as a linear system in ]K(m)[F]. This 

yields rational functions hj{u) G K{u) such that, for / € {1, . . . ,%}, yi — hi{u) = 0. 
7: Return the rational parametrization 

g{u)=0 

Xi=g\{u) y^=hi{u) 

Xn,-l = gn,-l{u) yn,-\ = Ky-\{u) 

Xn, = M - ^ CLlgl{u) yn, = hn,{u) 
f=l 



We now prove that the complexity of Algorithm [T] is bounded by the complexity of the underlying 
generalized MinRank problem and that most choices of (tti, . . . , o:„^-i) and M do not fail. 
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Theorem 6. Let /i, . . . ,fn,+ny S ]K[X,F] be an affine system of bi-degree (D, 1) such that the ideal 
(/i , . . . ,fny+ny) radical and 0-dimensional. Then there exists non-identically null polynomials h\ G 
. . . and h2 G K[zi, . . . ,Zny{n,.+n^.)] ^^ch that, for any choice o/ (tti, . . . ,a„^._i) and M = 

(mij) e ]K">' verifying: 

• the ideal (/i , . . . ,fn,+ny) i^ proper, radical and 0-dimensional; 

• the matrix jacy (/i , . . . , /n^+n,, ) verifies the conditions of Theorem 

• /Ji («!,... , a„^_i)/i2(mi J, . . . ,m„^,_„^+„^,) / 0, 

Algorithm^retums a rational parametrization of the variety of the system and its complexity is upper 
bounded by 

o + (^("^ + + y + (d- + ' 

Proof. Let / denote the ideal generated by /i, . . .,fn^+ny- According to ll22l |2l, for any radical 0- 
dimensional ideal, there exists a polynomial h^ such that if h\ {a\ , . . . , «„^-i ) 7^ 0, then the system is 
in shape position after the change of coordinates 

The polynomial h2 is chosen such that if hzimij) / 0, then the linear system f\ =•••=/«,.= 
in lC(M)[y] has rank exactly ny. Therefore h2 can be chosen as the GCD of the coefficients in 
the ring . . . ,z„^.(„^_,_„^)] of all terms u^ in the determinant of the following linear system in 

K[zi,i, . . . ,Zny,n,+ny , "] [Y] (whcrc the variables are 3^1 , . . . ,3;„p: 



^Zl,l ••■ Zl,n,.+nv 



/ fiig\{u),.-.,gn,-i{u),u,yi,...,yn,) modg(M) 



\Zny,l ■•• Zn,,n,+",/ (§1 1 ("), ",^1 jjn,) mod ^(m) 

Therefore this determinant is not null if and only if h2 does not vanish. Since the ideal generated by 
the input system (/i , . . . , ) is 0-dimensional and proper, it has finitely-many solutions, and hence 
/j2 is not identically zero. 

• The complexity of the substitution to compute the polynomials /, is upper bounded by 0{(nx + 
ny)Dnxny). 

• By Theorem [5| the complexity of the Grobner basis computation is upper bounded by 

• Since deg(^„J < DEG(/), a monomial u"'Y\'i^i^ xf' of de gree D can be evaluated in the uni- 
variate polynomials {gi{u),. . . ,gn,-i{u)) modulo g{u) in complexity C?(DDEG(/)) by using a 
subproduct tree |4|, quasi-linear multiplication of univariate polynomials and quasi-linear mod- 
ular reduction. Since there are at most {n^ + + 1) ("'^^) such monomials in the system 

/i , . . . , /„,+„,, , the Step 4 of the Algorithm needs at most O (n^ + ny)ny ^) D D EG (7)^ arith- 
metic operations in K. 

Notice that + < (^^l) and DEG(/) < (^^"'t/'^')- 
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- If D > 2: for any a,b,c en such that b < a, Qc < Therefore, DuyC'^f) < 
^n,-+ny+2Dy ^j^q^ notice that, for D > 2 and for any nx,ny such that rixny > 1, fix + riy + 
2D < D{nx + fly) + 1. Therefore, 

- IfD= 1: + = (n, + is bounded by i:'i;:"0 i^'^';^^') ■ 

Therefore, the complexity of the Step 4 of Algorithm [Tj^is upper bounded by the complexity of 
the Grobner basis computation: O ) • 

• To solve the linear system by using Cramer's rule, we need to compute + 1 determinants of 
{rix X «v)-inatrices whose entries are univariate polynomials of degree D. This can be achieved 
by using a fast evaluation-interpolation strategy with complexity C?(D«f+') (since multi-set 
evaluation and interpolation of univariate polynomials can be done in quasi-linear time, see e.g. 
11). 

Since DEG(/) is bounded by D'^'^ {"'n"') > these complexities is upper bounded by 

H(:-T)r'::'")"--(-'("':"'))> 

□ 

Remark 4. According to / TO] Lemma 15] and / TO] Lemma 16], if D = 1, there exists a non-empty 
Zariski open subset 0\ of the set of systems of bi-degree (1,1) which are 0-dimensional and radical. 
The proofs are similar for systems of bi-degree (D, 1) with D G N. 
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